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Abstract 

We deal with the elastic scattering by a large number M of rigid bodies, Dm ^Bm + Zm, of 
arbitrary shapes with 0 < e << 1 and with constant Lame coefficients A and /r. We show that, when 
these rigid bodies are distributed arbitrarily (not necessarily periodically) in a bounded region Q, of R® 
where their number is M := M(e) := 0(e“^) and the minimum distance between them is d d{e) ~ e* 
with t in some appropriate range, as e —>■ 0, the generated far-held patterns approximate the far-held 
patterns generated by an equivalent medium given by uPplz — {K + l)Co where p is the density of the 
background medium (with I 3 as the unit matrix) and {K + l)Co is the shifting (and possibly variable) 
coefficient. This shifting coefficient is described by the two coefficients K and Co (which have supports 
in Q.) modeling the local distribution of the small bodies and their geometries, respectively. In particular, 
if the distributed bodies have a uniform spherical shape then the equivalent medium is isotropic while 
for general shapes it might be anisotropic (i.e. Co might be a matrix). In addition, if the background 
density p is variable in and p = 1 in R® \ II, then if we remove from Q, appropriately distributed small 
bodies then the equivalent medium will be equal to in R®, i.e. the obstacle characterized by p 
is approximately cloaked at the given and fixed frequency u>. Elastic wave scattering, Small-scatterers, 
Effective medium. 

2000 Math Subject Classification: 35J08, 35Q61, 45Q05 


1 Introduction and statement of the results 
1.1 The background 

Let i?i, ^ 2 , ■ •., Bm be M open, bounded and simply connected sets in with Lipschitz boundaries, con¬ 
taining the origin. We assume that their sizes and Lipschitz constants are uniformly bounded. We set 
Dm '■= eSm + Zm to be the small bodies characterized by the parameter e > 0 and the locations 
TO = 1,..., M. 

Assume that the Lame coefficients A and p, are constants satisfying p > 0 and 3A -I- 2p > 0 and the 
mass density p to be a constant that we normalize to a unity. Let U’’ be a solution of the Navier equation 
(A® -I- u}'^)W = 0 in A'^ := {pA -|- (A -b p)V div). We denote by C/® the elastic field scattered by the M 

small bodies C due to the incident field C/L We restrict ourselves to the scattering by rigid bodies. 

Hence the total field [/* := W + satisfies the following exterior Dirichlet problem of the elastic waves 


(A^-bwAV* =0 inKA ( U Dm], 

(1.1) 

\rn—l J 

U*\3D^ = 0, 1 < TO < M 

(1.2) 
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with the Kupradze radiation conditions (K.R.C) 


lim 

|t|—)- oo 


, dUp 

'5|a;| 


itvpUp) — 0 , 


and 


lim 

|a;|—>-oo 



IKsUs) = 0 , 


(1.3) 


where the two limits are uniform in all the directions x := S S^. Also, we denote Up := ^V(V • U^) to 
be the longitudinal (or the pressure or P) part of the field and Us '■= x (V x C/®) to be the transversal 
(or the shear or S) part of the field corresponding to the Helmholtz decomposition = Up + Us- The 
constants Kp and Ks '■= ^ are known as the longitudinal and transversal wavenumbers, Cp := \/A + 2/r 

and Cs := y/Jl are the corresponding phase velocities, respectively and w is the frequency. 

The scattering problem ( l.lfO ) is well posed in the Holder or Sobolev spaces, see [Colton fc Kress, 1998| 
Colton fc Kress, 1983| Kupradze, 1965 Kupradze et ah, 1979| for instance, and the scattered field U^ has 


the following asymptotic expansion: 


W{x) := 


olKpIxl 


piK,s\x\ 1 

+ O(^), |x| ^ 00 


(1.4) 


uniformly in all directions x £ S^. The longitudinal part of the far-held, i.e. 17“ (x) is normal to while 
the transversal part 17“(x) is tangential to S^. We set 17“ := (17“,17“). 

As usual, we use plane incident waves of the form W{x,9) := -I- giK,e-x^ where 9^ is any 

direction in perpendicular to the incident direction 9 G S^, a,/3 are arbitrary constants. The functions 
U^{x,9) := 17“ (x) and U^{x,9) := 17“ (x) for (x,0) € x are called the P-part and the S-part of the 
far-held pattern respectively. 


Definition 1.1. We define 

1. d := min dmj, wheredmj '■= dist{Dm, Dj). We assume that 0 < d < dmax, and dmax is given. 

2- Wniax as the upper bound of the used frequencies, i.e. oj G [0, Wmax]- 

5. n to be a bounded domain in containing the small bodies Dm, m = 1,..., M. 


1.2 The results for a homogeneous elastic background 

We assume that Dm = cBm + Zm, m = 1,..., M, with the maximal diameter a := emax^ diam(Hm), are 
non-hat Lipschitz obstacles, i.e. are Lipschitz obstacles and there exist constants to G (0,1] such 

that B^^a{zm) C Dm C B\{zm), where tm are assumed to be uniformly bounded from below by a positive 
constant. In jChalla &: Sini, 2015| , we have shown that there exist two positive constants eg and Cq depending 
only on the size of D, the Lipschitz character of Bm, to = 1,..., M, dmax and oJmax such that if 

e < Co and VAf — 1- < cq (1-5) 

then we have the following asymptotic expansion for the P-part, 17“(x, 0), and the S-part, 17“(x,0), of the 
far-held pattern: 


C/“(x,d) = 


+0 M 


iTTCp 


(x 0 x) 


■ M 

E 

m—l 


d5-3a d9-6a 


-l — X-Zm. y-\ 

e Qn 


M(M- 1) 




+ M(M-1)2^ 


( 1 . 6 ) 
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U^{x,9)= - —^(I-i(g)i) 

Air ci 


M 

E 


— I - X- Z-n 


Qn 


.m—1 


+0 M 


d5-3a 


M{M -1) 


g3 ^4 


(P°‘ 


M{M -if 


<P° 


(1.7) 


where a, 0 < a < 1, is a parameter describing the relative distribution of the small bodies. 

The vector coefficients Qm, m = l,M, are the solutions of the following linear algebraic system 


M 


C-^Qm = 


( 1 . 8 ) 


i=i 

j^m 


for m = 1,..., M, with r‘^ denoting the Kupradze matrix of the fundamental solution to the Navier equation 
with frequency w, Cm ■= Jgj^ <Xm{s)ds and am is the solution matrix of the integral equation of the first 
kind 


I r (^Sm ^ ^}^m(^)ds — I, Sni € ^Dn 

JdD^ 


(1.9) 


with I the identity matrix of order 3. 


Consider now the special case dminf* ^ d < dmax^* and M < Mmax^ ^ with t, s > 0, dmin, dmax and 
Mmax are positive. Then the asymptotic expansions (1.6|1.7) can be rewritten as 


M 




-l — X-Zjj 


Qr 


m—1 


s _|_ ^3—s — 5t-\-3ta _|_ ^4—s—9i+6ia _|_ ^3—2s—2ta _|_ ^4—3s — 3ta _|_ ^4—2s—5i+2ia^ j 
M 


1 




^ m—1 

s _|_ ^3—s — 5t-\-3tot _|_ ^4—s—9i+6ia _|_ ^3—2s—2£a _|_ ^4—3s —3£a _|_ ^4—2s—5£+2£a^ j 

As e —> 0, the error term tends to zero for t and s such that 

5 3 4 

0 < t < 1 and 0 < s < min{2(l — t), 3 — 5< + 3ta, 2 — -t + ta, 4 — 9t + 6ta, - — ta, - — to}. (1-12) 
In [Challa fc Sini, 2015|, we have shown that Qm ~ e, then we have the upper bound 


M 


I V e-^---^-Qm\ < M sup \Qm\ = 0(e'-9- 

1 m—l,...,M 

m—1 ' ’ 


(1.13) 


Hence if the number of obstacles is M := M(e) := 0(e ®), s < 1 and t satisfies (1.12), e —>■ 0, then from 

(1.14) 


(1.10 1.11), we deduce that 

U°°{x, 0) —>■ 0, as e —0, uniformly in terms of 9 and x in S^. 


This means that this collection of obstacles has no effect on the homogeneous medium as e —>■ 0. 


Let us consider the case when s = 1. We set H to be a bounded domain, say of unit volume, containing 
the obstacles Dm,rn = 1,...,M. Given a positive and continuous function AT : —)■ K, we divide H into 
[e“^] subdomains flm, m = 1,..., [e“^], each of volume with Zm G as its center and contains 

[K{zm) + 1] obstacles, see Fig[^ We set Kmax ■= sup^^{K{zm) + 1), hence M = Y}f^i[K{zm) + 1] < 
Kmaxf~f = 0(e“^). 
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Figure 1: A Schematic example on how the obstacles are distributed in fl. 


Theorem 1.2. Let the small obstacles be distributed in a bounded domain Ll, say of unit volume, with 
their number M := M(e) := 0(e“^) and their minimum distance d := d{e) 0(e*), ^ < t < ^, as 

e —>■ 0, as described above. In addition, we assume that the Lame coefficients A and ^ satisfy the conditions 
max{«„ Kp} < and c(A + 2/r) < max^c-(a„) Q 

1. If the obstacles are distributed arbitrarily in LI, i.e. with different capacitanees, then there exists a 
potential Cq G np>iL^(K^) with support in LI such that 

lim U°°{x,0) = Ulff{x,9) uniformly in terms of 9 and x in (1-15) 

a—>-0 

where {x, 9) is the farfield corresponding to the scattering problem 

(A^+a;^-(h:+l)Co)17,5 = 0 (1.16) 


UoldDrr, = 0 , 1 < m < M 


with the radiation conditions 


lim 

|ai|—>-oo 


kK 


dU, 


0 ,p 


cl|a:| 


'iripUo^p'j — 0 , 


and 


lim 

|ai|—>-oo 


kl( 


dU, 


0,s 


a|x| 


iKsUo,s) = 0 


(1.17) 


(1.18) 


2. If in addition K |a is in C'°’'>'(f2), 7 S (0,1] and the obstacles have the same capacitances C, then 

U°°{x,9) = U^{x,9) + 0 (e™“{'>'’ 3 ’i“ 3 t}^ uniformly in terms of 9 and x in (1-19) 

where Cq = C in LI and Cq = 0 inM.^ \ LI. 

^The constant c is defined in and denotes the capacitance of each scatterer (i.e. the acoustic capacitance). These 

conditions on the Lame coefficients A and p. can be replaced by a condition on the wave number u, see Lemma |2.3| and the 
related footnote. 
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1.3 The results for variable background elastic mass density 


Assume that the Lame coefficients A and fjL are constants satisfying fj, > 0 and 3A + 2/j, > 0 and the mass 
density p to be a measurable and bounded function which is equal to a constant that we normalize to a unity 
outside of a bounded domain fl. We set Pmax to be the upper bound of p. 

In this case, the total field C/* ■,= U'^ + Up satisfies the following exterior Dirichlet problem of the elastic 
waves 

Dr^ , (1.20) 

Ul\ 9 D^ = 0, 1 < m < M (1.21) 

with the Kupradze radiation conditions (K.R.C) 


(A^ + w^)!/* = 0 in 


M 

u 

m—1 


lim 

|ai|—>-oo 


,dU, 


p,p 


d\i 


iKpUp^p) = 0, and 


lim 

|ai|—>-oo 




IKsUp^s) = 0, 


( 1 . 22 ) 


where the two limits are uniform in all the directions * ]yy £ and Up^p and Up^s are respectively the 
P-part and S-part of the scattered fie ld ^ 

The scattering problem ( 1.20||1.22 ) is well posed in the Holder or Sobolev spaces, see 
[Colton fc Kress, 1998[ [Colton fc Kress, 1983| [Kupradze, 1965[ [Kupradze et ah, 1979[ for instance, and the 
scattered field t/® has the following asymptotic expansion: 


„iKp\x\ piK,s\x\ 1 

u;{x) := (i) + _[/- (i) + 0 (_), |a;| ^ OO (1.23) 

uniformly in all directions x S S^. The longitudinal part of the far-held, i.e. U^p{x) is normal to while 
the transversal part U^{x) is tangential to S^. We set := (C/^, U^^). 

As in the case of constant background mass density, there exist two positive constants cq and cq depending 
only on the size of H, the Lipschitz character of Bm, to = 1,..., M, dmax, Wmax, Pmax and Pmax such that if 

e < Co and \/M — 1- < cq (1-24) 

a 

then we have the following asymptotic expansion for the P-part, U^p{x, 9), and the S-part, U^g{x, 9), of the 
far-held pattern: 


U^p{x,9)= K~(i,0) + 

+o(m 


M 


G'^p{x,Zm)Qp,r. 


.m—l 


^3 ^4 


(l5-3a 

M 




^ J5-2 






17“(x,0)= K“(x,0) + 

^3 

+0 ( M ' '' 


Y ‘^y!six,Zm)Qp,r. 


.m—l 
4 1 


^5-3a d^-Ga 


■M(M- 1) 


r e3 


A 1 


d2a (;5-2a 


■M{M -if 


d^° 


(1.25) 


(1.26) 


where G^p{x, Zm) and G^g{x, Zm) are the P-part and S-part of the farhelds of the Green’s function Gp(x, z), 
of the operator A® -|- uj^p in the whole space evaluated in the direction x and the source point Zm- 
The vector coefhcients Qp.m, to = 1,..., Af, are the solutions of the following linear algebraic system 


c-^Q 


p,m 


M 

^pi.Z'm’,9') ^ ^ G pi^Zm-, Zj)Q p.j ^ 

4 = 1 
j^m 


(1.27) 
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(1.28) 


for m = 1, with Vp{-,9) := Vp(-,6) + is the total field satisfying 

(A®+a;V)^p =0 in 

and the scattered field Vp{-,9) the Kupradze radiation conditions (K.R.C). 


Corollary 1.3. Let the small obstacles be distributed in a bounded domain 12, say of unit volume, with 
their number M := M(e) := 0{e~^) and their minimum distance d := d{e) 0(e‘), ^ < t < ^, as 

e —>■ 0, as described above. In addition, we assume that the Lame coefficients A and ^ satisfy the conditions 
max{K„ Kp} < and c(A + 2^l) < ■ f\ 

1. If the obstacles are distributed arbitrarily in LI, i.e. with different capacitances, then there exists a 
potential Cq S np>iL^(K^) with support in 12 such that 

lim U?°{x,9) = Uff^ffx,9) uniformly in terms of 9 and x in (1-29) 

where U^q{x, 9) is the farfield corresponding to the scattering problem 

(A^ + uj^p -{K + l)Co)17*_o = 0 m ]R^ (1.30) 

=0, l<m<M (1.31) 


with the radiation conditions. 

2. If in addition K |n is in C'°’'>'(r2), 7 G (0,1] and the obstacles have the same capacitances, then 

U^{x,9) = U^q{x,9) + 0 ( 6 ™™^'’'’uniformly in terms of 9 and x in (1-32) 

where Cq = C in LI and Cq = 0 m \ 12. 


1.4 Applications of the results and a comparison to the literature 

The main contribution of this work is to have shown that by removing from a bounded region of an elastic 
background, modeled by constant Lame coefficients A and p and a possibly variable density p, a number 
M := M{e) ~ of small and rigid bodies with diameter of order e, distant from each other of at least 
d := d{e) ~ e‘, g < t < A, then the ’perforated’ medium behaves, as e —>■ 0 , as a new elastic medium 
modeled by the same Lame coefficients A and p but with a new coefficient — {K + l)Co. The coefficient 
K models the local distribution (or the local number) of the bodies while the coefficient Cq, coming from 
the capacitance of the bodies, describes the geometry of the small bodies as well as their elastic directional 
diffusion properties (i.e. the anisotropy character). In addition, we provide explicit error estimates between 
the far-fields corresponding to the perforated medium and the equivalent one. From this result we can make 
the following conclusions: 


1. Assume that the removed bodies have spherical shapes. For these shapes the corresponding elastic 
capacitance C is of the form cl^ (i.e. a scalar multiplied by the identity matrix). In section 
we describe a more general set of shapes satisfying this property. Hence the equivalent coefficient 
uffp — {K + l)c is isotropic while for general shapes it might be anisotropic. To achieve anisotropic 
coefficients, a possible choice of the shapes might be an ellipse. 


^Same comments on c as in Theorem 


1.2 


apply here too. 
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2. If we choose the local number of bodies K large enough or the shapes of the reference bodies, 

m = 1, having a large capacitance (i.e. a relative large radius) so that uJ^p— {K + l)c < 0. This 

means that, at the frequencies ui for which this inequality is satisfied, the elastic material, located in 17, 
should have an unusual behavior, namely the traction forces (modeled by the Lame coefficients A and 
p) and the addition of the two body forces (modeled by the density p and the coefficient {K + l)Co 
respectively) will act in the same direction. 

3. Assume that the background medium is modeled by variable mass density p > 1 in 17. If we remove 
small bodies from 17 with appropriate K and/ or capacitance Cq so that p — {K + l)ca;“^ = 1, then 
the new elastic material will behave every where in as the background medium. Hence the new 
material will not scatter the sent incident waves at the given frequency w, i.e. the region 17 modeled 
by p will be cloaked at that frequency. 


The ’equivalent’ behavior between a collection of, appropriately dense, small holes and an extended 
penetrable obstacle modeled by an additive potential was already observed by Cioranescu and Murat 
ICioranescu fc Murat, 1982|[Cioranescu &: Murat, 1997] and also the references therein, where the coefficient 
K is reduced to zero since locally they have only one hole. Their analysis is based on the homogenization the¬ 
ory for which they assume that the obstacles are distributed periodically, see also [Bensoussan et ah, 1978 
and [Jikov et ah, 1994| . 


In the results presented here, we do not need such periodicity and no homogenization is used. Instead, 
the analysis is based on the invertibility properties of the algebraic system (1.81 and the precise treatment of 
the summation in the dominant terms of (1.6)-(1.7). This analysis was already tested for the acoustic model 
in [Ahmad et ah, 2015| . Compared to [Ahmad et ah, 2015| , here, in addition to the difficulties coming from 
the vector character of the Lame system, we improved the order of the error estimate, i.e. 
instead of 0(e™'"{'>'’3“l‘l) which, for t := | for instance, reduce to 0(e™“ffi’sl) and w}) respec¬ 

tively. 

Let us mention that a result similar to (1.29), for the acoustic model, is also derived by Ramm in several 
of his papers, see for instance [Ramm, 2007| , but without error estimates. Compared to his results, and as 
we said earlier in addition to the vector character of Lame model, we provide the approximation by improved 
explicit error estimates without any other assumptions while, as shown in [Ramm, 2007[ [Ramm, 201 1| for 
instance, in addition to some formal arguments, he needs extra assumptions on the distribution of the 
obstacles. 

We discussed above the application of our result in elastic cloaking via perforation (or perturbation with 
small bodies). Another way of achieving elastic cloaking is by using transformations, see [Hu &: Liu, 2015| . 

The rest of the paper is organized as follows. In section 2, we give the detailed proof of Theorem [1.2| 
In section 3, we describe the one of Corollary [1.3[ by discussing the main changes one needs to make in the 
proof of Theorem 1.2 Finally, in section 4, we discuss some invariant properties of the elastic capacitance 


to characterize the shapes that have a ’scalar’ capacitance. 


2 Proof of Theorem 1.2 


2.1 The fundamental solntion 

The Kupradze matrix F"^ = (r/))/^^;^ of the fundamental solution to the Navier equation is given by 
r‘^{x,y) = -$«^(a;,2/)I-h ^VxVj[l>„^(x,y) - $„p(x,y)]. 




( 2 . 1 ) 


.iK\x-y\ 


where $^( 0 ;, y) = denotes the free space fundamental solution of the Helmholtz equation (A-I-k^) u = 

0 in K^. The asymptotic behavior of Kupradze tensor at infinity is given as follows 


F-(x,j/) = 


^ giKpl^l 

X ® X —;—; 1 


iTTCp 


47r ci 


ski 


■y + o{\x\-^) 


( 2 . 2 ) 


7 







































with ^ = if]- G and I being the identity matrix in see Alves fc Kress, 2002| for instance. As mentioned 
in Ammari et ah, 2007|, (2.11 can also be represented as 


- OO ./ ^ 

1 1 




1 A (;-i) 


drr l\[l + 2) 

from which we can get the gradient 




) \x-y\'‘ ^{x-y)®{x-y), 


(2.3) 


[l-l) 


V.r-(x,y) = E + 2) ~ ~ 


3 3 


- (4^^k-l/l' ^ I (/-3)|a;-?/| 2(g)3(a._y)+i0(x-y) + ^^(x-y)jej/(g)ej®e 

j=i j'=i 


where Cj is standard basis vector. Using the formulas (2.3) and (2.4) we can have the following estimates, 
for x,y G ^},x y, see [Challa fc Sini, 2015]; 


with 


|r“(x,y)|<^ 


Cj 


■-y\ 


■Ch 


|v,r-(x,2 /)|< — 


Co 


'- 2 /p 


C 
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(2.5) 


^7:= 


Cg := 2^ 


a; / 1 — i^Ks diam(U)) 


1 1 

72 72 

Nn 


1 — diam(f2) 2^^ ^ 


\ uj f 1 — (j/Cp diam(fl))^“^ 
/ c3 I 1—iKpdiam(U) 


1 


/ 1 1 — (|ks diam(r2)) 


Nn 


^ rf 24 


1 — ^Kg diam(fl) 2^^‘ ^ 


/ 1 1 — diam(il)) 

5—r ' 


2Nn-l 

Nn 


20 


1 — ^Kp diam(U) 2^^ ^ j ’ 


and Nq = [2 diam(U) maxims, Kp}e^] where we assume the frequency ui and the Lame parameters A and /i 
to satisfy the condition max{Ks, Kp} < . 

Indeed, for x G Dm and s G Dj, first let us recall the following result from |Challa & Sini, 2015[ Lemma 

2 . 6 ]. 

Lemma 2.1. For each k > 0 and for every n G Z"*" with n > ke^ [=: /V(fc)] we have n\ > 


Now, by making use of (2.3), we get the following estimate; 
1 1 


|r^(3^,2/)l < 


{kI + kD \x-y\~ 


Fk uj‘ 

_ 1 

A'TT a j ( 1 — 


1 


( 2 , 


^/+2 I j ^ Z +2 


< - 


47r (I — !)!(/ + 2) 

—I (ks + «p) + ^ 7771 


1 1 
47r uF 


p )\x-y\ 
1 


Z-1 


- 2/1 


1^1 


[I-mi+ 2 ) 




^-1 










































< — 


1 

47r 


- y\ 


„2 + „2 ) + X! - 


^ G-l)!a + 2 ) 




< 


[By writing Nq = [2 diam(r2) max{Ks, Kp}e^] and using Lemma |2.l| 


1 

4:TT 


- y\ 




LU 


/ Nfi X - \ t — ± -1 


> \i=i 


2/-1 

l=Nn+l 
1—1 00 


l=Nn+l 


1 

47r 


1 

4:TT 


\x-y\ Ws cl 


,Nn 


Ks I I - (|Ksdiam(Sl)) “ 1 


' 1— iKsdiam(r2) ' 2^f* ^ 


Kp (l - (dKpdiam(ri))^‘^ 1 

Cp I 1 — jKp diam(ri) 


C 7 


- y\ 


+ Cs 


Now, by making use of \2A^ we get the following estimate; 

|V,r‘"{a:, 2 /)| < 


1 1 

477 


1 


7 {kI + Kl)\x-y\ 2 + - (14^4 + 12k^) 


1 


00 

-—V- 

47r^Z(Z-2)!(; + 2 ) 


((2Z + 8)Ks"''^ + (/ + 8)7vp''"^) [a; — y\^ 


-2 


1 1 

~ 477 07^ 


< - 


. _ y \2 + '^P) 


+E 


pt ■ 4 

1 


Z=2 


Z(;-2)!(Z + 2) 


(( 2 ; + 8 )k£^ + (? + 8 )Kp+^) diam(f 2 )^ 


1 



1 ^ 

+ 1 

^07^ a; 2 \ 

477 

Jx - y |2 


' cIj 

+ 4 

V cf 4 ) 


00 1/2 2 \ 

+ E ^i_ 2)1(1+ 2) ^ " 


[By recalling Nq = [2 diam(f2) max{Ks, Kp}e^] and using Lemma 


2.1 


< 


1 

477 


7 


-y|^ 


1 1 

^ 72 

‘-p 


1 . o;^ 

Cp 


+6^ f E fxK«diam(L!) ) + E (b 


z=o 
2 / ^£1 — 1 


I 00 / ^ \ /^ 


p \ z=o 


+®74 E x«pdiam(n) + U 


l^Nn 
I 00 



l=Nn 



Co 


477 [ [x — 7/1' 


+ C: 
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The estimates in (2.5) can be written as 

|r‘^(^, 2 /)l< " 


47r| X — 7/1 ’ 


|Vpr-(x,7/)|< 


47r| X — 7/|2 ’ 


( 2 . 6 ) 
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for different points x,y ^ VL, where 


c := maxICr, Cg diam(ri), Cg, Cio diam(r 2 )}. 


(2.7) 


2.2 The relative distribution of the small bodies 

The following observation will be useful for the proof of Theorem |1.2[ For m = 1,..., M fixed, we distinguish 
between the obstacles Dj, j ^ to by keeping them into different layers based on their distance from Dm- Let 
us first assume that K(zm) = 0 for every Zm- Hence each Dm bas the (same) volume e and contains only one 
obstacle Dm- Without loss of generality, we can take the H^’s as cubes. Hence we can suppose that these 
cubes are arranged in cuboids, for example Rubik’s cube, in different layers and Dm is being located at the 
center. Observe that the total number of cubes in a Rubik’s cube consisting of N layers is (27V + 1)^ cubes, 
see Fig[^ Since, the volume of each Dm is e and the cubes are arranged in N layered Rubik’s cube, the total 


number of cubes are e ^ which is then equal to (2A^ +1)^ and then N = 


' 3-1 

2 


. Hence, (1) the total cubes 


upto the layer consists of ( 2 n + 1 )^ cubes and ( 2 ) the number of obstacles located in the n 7 ^ 0 , layer 

/ ^ max diam(B.^)'' 

is [(2n + 1)^ — {2n — 1)^] = 24n^ + 2, and their distance from Dm is more than n I es — e- 


■3-1 


for n = 0 ,..., 

Now, we come back to the case where K{zm) ^ 0. First observe that \ < 1. Hence with 

such Hm’s, the total cubes located in the layer consists of at most the double of [( 2 n + 1 )^ — ( 2 n — 1 )^], 
i.e. 48n^ + 4. It is due to the fact that upperbound of the fraction relates to the lower bound of 

the total number of cubes ( which is explained in the previous paragraph for the case fl^’s of volume e) and 
the lower bound of the fraction relates to the upper bound of the total number of cubes located 

in each layer, which is 2 [( 2 n + 1)^ — ( 2 n — 1)^]. 



Figure 2: Rubik’s cube consisting of a single layer 


2.3 Solvability of the linear-algebraic system (1.8) 


We start with the following lemma on the uniform bounds of the elastic capacitances, see 
[Challa fc Sini, 2015[ Lemma 3.1] and the references therein. 

^Indeed, if we take any two adjacent cubes and f22 which have the obstacles Di and D 2 at their centers respectively, then 


the distance between the obstacles is more than 
and the volume of the cube Qm is e. 


max diam(Bm,) 
l<m<M 


as the maximal radius of Dm is e ^^ 


max diamlBrn,) 
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Lemma 2.2. Let A™" and A™“^ be the minimal and maximal eigenvalues of the elastic capacitance matrices 
Cm, for m = 1,2,... ,M. Denote by the capacitance of each scatterer in the acoustic casej^ then we 
have the following estimate; 

piC'^{Bm)e = Pi AS“ < XTgZ ^ Cf^ = {X + 2p)C^{Bm)e, (2.8) 

for 771 = 1,2,..., M. 


The constant C{Bm) is the acoustic capacitance of the reference body Bm which can be estimated above 
and below by the Lipschitz character of Bm, [Challa fc Sini, 2014|. T he following lemma provides us with 
the needed estimate on the invertibility of the algebraic system (1.8) whose coefficient matrix ’B’ is given 
by; 

/ -r‘"(0i,z2) -r‘"(zi,z3) ••• -r-(zi,0M) \ 

-r“(z 2 , 2 i) -r“(z2,23) ••• -r-(z 2 ,^M) 


B := 


\ -T'^{zm,Zi) -T^{zm,Z 2 ) 


-T^{zm,zm-i) 


-c 


(2.9) 


M 


/ 


Lemma 2.3. The matrix B is invertible and the solution vector Qm rn = 1,...,M, of (1.8) satisfies the 
estimate: 

M 

( 2 . 10 ) 

\ TT / --- 

m—1 


E ll<3,»lli < ((A + 2 rt-‘ - ^ 

m—1 \ / m—1 

ifmaxCZ e-^ < <^ddztion, 


M 


Y. WQrnWl < ( (A + 2p)-^ - c(max^C°) 


-1 


n—1 


maxC^ M max||[/*(2:m)||2. (2.11) 

m—1 


0 


Since Cm = C°'{Bm)c, then 


makes sense if m.aXmC°‘{Bm)M„ 


maxCf 


e-i < 


Vl^Mmax c(A + 2 /i) 


( 2 . 12 ) 


< 


^T4c(A+2/j) 


. As C°‘{Bm) is proportional to the radius of Bm, then 


(2.12) will be satisfied if uj and the Lame parameters A and p, satisfy the condition 


maxji 


J < 


diam(n) 


and c(A + 2p) < 


VliMmax maXmC'“(B„) 


(2.13) 


■^Recall that, for m = 1,. .. ,M, := Jq^ (Jm{s)ds and <Jm is the solution of the integral equation of the first kind 

Jqd 4^'|t-sl = 1) ^ £ dDm, see [Challa & Sini, 2014| . 

^ Alternatively, see [Challa fc: Sini, 2015| Lemma 3.2], we can have the following estimate (and then the invertibility of B): 


3t 


XlllQm ||2 <4 ( -3^ max A) 


mill -max p \ max ) "N 

l<m<M 57rdl<m<M J 

^ ' m = l 


\4 M 




l<m<M 


if ( max )<r min [with the positively assumed value 

eia^ I \ 3 l<m<M ^^9^1 


^—2diam{Q) - 


1 — (^ K. s d i a m (f 2 


diamm) +^Nk^]-diam{n) j 


1 — (^ K,pdiam(Q,'j^ 

1 — (^ Kpdiam{Q,'jj 



with max{/ts, K,p} < Observe that this last condition and the positivity of r can be seen as conditions on the wave 

number uj. Using these conditions, we can replace in Theorem [1.2[ and Corollary |1.3[ the conditions on the Lame coefficients A 
and p. by a condition on u. 
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recalling that c is defined in (2.7). Finally, let us observe that the right hand side of (2.13) depends only 
Mmax and the Lips chit z character of the reference obstacles Em’s. 

Proof of Lemma \2^ We start by factorizing B as B = — (C^^+B„) where C := Diag{Ci,C 2 , ■ ■ ■, Cm) G 
]^Mx M^ j jg identity matrix and B„ := —C~^ — B. We have B : ^ so it is enough to prove 

the injectivity in order to prove its invertibility. For this purpose, let X, Y are vectors in and consider 
the system 


(c-i + B„)x = y. 


(2.14) 


We denote by and (.)®'"® the real and the imaginary parts of the corresponding complex vector/matrix, 

we also set by C/. From (2.14) we derive the following two identities: 

(C7+B;“')X’'“'-B™9y™9 = ^2.15) 

(C 7 + B;“')X*'"s + B™9X’'“' = (2.16) 


and then 


^ ■Qreal'^j^real j^real'j ^-^irngj^img j^real'^ _ ^ y^reaZ j^real'j 

^ j^img\^ _j_ ^ -Qimg ^real j^img'j _ ^ yrimg j^img'j 


Summing up (2.17) and (2.18) we obtain 

’'real vreal\ 


^reaL\^ _|_ ^ greaZj^rea/ j^real^ _|_ ^ qa^^ img ^img\^ _|_ ^ -Qreal^img j^img'^ 

_ l^^real ^real\^ _j_ ^ ^img\^ 


(2.17) 

(2.18) 


(2.19) 


( 2 . 20 ) 


The right-hand side in (2.19) can be estimated as 

^ j^real j^real^^l/2 ^y^real yrreal^^l/2 ^ j^img j^img^^l/2 ^y^img Yimg'<^lj2 

< 2(yl'^xl'l)i/2(rl'l^yl'l)i/2_ 

Here, := (|jX’'®“^|p -|- ||X*'"®|p) 2 . Let us now consider the right hand side in (2.19). First we have 

|( ( 2 . 21 ) 

where ||H™'||i := and := r“(z„z,) if z ^ j and (B^^Om - « io^J = 

1,...,M. Hence from (2.3) ; * 7 ^ j- From the observation before Lemma 

deduce that 


2.2 


we 


M 


E iBr% < mm„ 

i ,1 = 1 


< MM„ 


c 2 


(47r)2d2 


-f J 2 2 [( 2 n-f l)3-(2n-l)3]- 


(47r)2n2 


(47r) 


_ 2 |e 3 
.9 € 3 

13c2^ 

TT^ 

n—1 


El 


+ e' 

t2 


( 2 . 22 ) 


or 


and then 


\-nreal\\ ^ l^max C 
2 < - 


(13e 


-2 , 


(2.23) 


real \rreal -\rreal\ 


Qreai 


^rea()| |( X*”*®) | < 




-(13e 


-2 I ,- 2 i-l^i 


1^1 


(2.24) 
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Using Lemma 2.2 we deduce that 


C^J^rea^^^rea/) (^Q i > ( A + 2^)" ^ (maxC^) |X 


From (2.19), (2.20), (2.24) and (2.25), we deduce that 


M 


-2 


M 


E < ( (A + 2M)--(maxC;.)-- - |y„. 


(2.25) 


(2.26) 


if (A + 2/.t) ^{maxC^) ^ > 
to be 


which holds true only if t < A. Hence, (2.261 turns out 


M 


-2 


M 


E < ( (^ + 2M)-'(maxC“ )-i - 


(2.27) 


if (A + 2/r) ^(maxC)),) 


-1 ^ %/l4Mmax C£~^ 


and then the matrix B in algebraic system (1.8) is invertible. 


2.4 The limiting model 

From the function K, we define a bounded function Kf : 


as follows: 


K fr) ■= K (z ) ■= I + 1 if X € flm 

e\ ) ■ e\ m) ■ Q jf j. ^ fgj. g^jjy rn = I,. . . , [a“^]. 


(2.28) 


Hence each Hm contains [K^{zm)] obstacles and K^ax '■= sup^^ K^{zm)- 

Let Ce be the 3x3 matrix having entries as piecewise constant functions such that Ce|n,„ = Cm for all 
m = 1,... ,M and vanishes outside fl. Here, Cm are the capacitances of B^’s. From [Challa fc Sini, 2015| , 
we can observe that Cm are defined through defined through Cm ■= Cm^, and are independent of e. 

We set 

(2.29) 


C := max HCmlloo- 

l<m<M 


Consider the Lippmann-Schwinger equation 

Ye{z)+ [ T‘^{z,y)K^{y)C^(y)Y^(y)dy = U\z,0),zGn 

Jn 

and set the Lame potential 

V{Y){x):= [ T^{x,y)K^{y)C,{y)Y{y)dy, xGM.^. 

Jn 


(2.30) 


(2.31) 


The coefficients iFe and are uniformly bounded. The next lemma concerns the mapping properties of 
the Lame potential. These properties are proved for the scalar Poisson potential in [Colton fc Kress, 1998| , 
for instance. Similar arguments are applicable for the Lame potential as well, so we omit to give the details. 

Lemma 2.4. The operator V : L'^{n) -G is well defined and it is a bounded operator for any bounded 

domain f2 inM.^, i.e. there exists a positive constant Cq such that 


l|T(T)||ii-2(n) < co||T||L2(n). 


(2.32) 


We have also the following lemma. 
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Lemma 2.5. There exists one and only one solution Y of the Lippmann-Schwinger equation [ 2.30) and it 
satisfies the estimate 


|>"||L-(a) < CWWhhu) and ||Vr||i»(n) < C'W\\ 


(2.33) 


where being a large bounded domain which contains fl. 
Proof of Lemma 1^.51 

we see that I + V : Lf{Ll) 


Using the Lemma 


2.4 


can apply the Fredholm alternative to / + U : Lf(Ll) 


L^{n) is Fredholm with index zero and then we 
L^(n). The uniqueness is a consequence of the 


uniqueness of the scattering problem corresponding to the model 


(A®= 0, in (2.34) 

where Y := U* + Y^ and F® satisfies the Kupradze radiation conditions and F* is an incident field. 

The estimate ( |2.33 ) can be derived, as it is done in [Ahmad et ah, 2015| for the acoustic case, by coupling 
the invertibility of / + F : L'^(Lt) i—>■ L‘^(Lt) and the interior estimates of the solutions of the system 

(A'=+a;2/-iF,C,)F = 0. 

□ 


2.4.1 Case when the obstacles are arbitrarily distributed 

The capacitances of the obstacles Bj, i.e. Cj are bounded by their Lipschitz constants, see 
jChalla fc Sini, 2014| , and we assumed that these Lipschitz constants are uniformly bounded. Hence 
is bounded in L^(H) and then there exists a function Cq in L^(H) (actually in every L^’(H)) such that 
converges weakly to Cq in L^(H). Now, since K is continuous hence ATe converges to {K + 1) in and 

hence in Lf{Ll). Then we can show that converges to {K + l)Co in Lf{Ll). 


Since iFCe is bounded in then from the invertibility of the Lippmann-Schwinger equation and 

the mapping properties of the Lame potential, see Lemma 


2.5 


we deduce that |jl7*||i/2(Q) is bounded and 
in particular, up to a sub-sequence, C/* (:= F^) tends to Uq in L^(f2). From the convergence of to 


{K + l)Co and the one of U* to Uq and (2.30), we derive the following equation satisfied by Uq{x): 


U*{x)+ [ T^{x,y){K + l)CQiy)U*{y)dy = W{x,e) in H. 
Jn 


This is the Lippmann-Schwinger equation corresponding to the scattering problem 

(A'= + uj'^-{K + l)Co)U*Q = 0 in 


(2.35) 


(2.36) 


with Uq = Uq + W , and U® satisfies the Kupradze radiation conditions. As the corresponding farfields are 
of the form 

= - [ -^{x^x)e-^^^-^K+l)CQiy)U*Qiy)dy, 

Jn Cp 

USffix,0) = - f ^(I-i®i)e-*^^'^(iF+l)Co(j/)C/‘(r/)dy 

Jn Cg 

and the ones of U^ are of the form 

Uffp{x,9) = - [ -^{x(^x)e-^^^-^K,CMUKy)dy 

Jn 4’’" Cp 

C/“(x,0) = - [ ^(I-x^x)e-^^^-yK,CMUKy)dy 
Jn 4’’" Cs 

we deduce that 


Ufff{x, 9) — U^p{x, 9) = o(l) and Uff.{x, 9) — U^g{x, 9) = o(l), e —0, uniformly in terms oi x,9 S 
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2.4.2 Case when K is Holder continuous 

If we assume that K S 7 S (0,1], then we have the estimate \\{K + 1) — , e « 1. 

Since the capacitances of the obstacles are assumed to be equal, we set Cg to be a constant in fl and Cg = 0 
in \ n. Recall that C/g and are solutions of the Lippmann-Schwinger equations 

Uo + f (K + l)T‘^ix,y)Coiy)U^oiy)dy = 

Jn 


and 

U,+ f K,T^ix,y)Co{y)Ul{y)dy = U\ 

Jn 

From the estimate \\{K + 1) — < Ce^, e « 1, we derive the estimate 

U§°{x,9) — U^{x,6) = 0{e^), e « 1, uniformly in terms of i, d € S^. 

2.5 The approximation by the algebraic system 

For each m = 1,..., M, we rewrite the equation ( |2.30[ ) as follows 

M 

UejZm) + (.^mi Zj)CjUe{Zj)e = U’'{Zm, 9) + A + B , 


i=i 

j^m 


where 


’ - r 

A:= ^r‘^(zm,Zj)J^e(zj)CjU^(zj)Vol(nj) - / T‘^{zm,y)K^{y)C^{y)U^{y)dy 

.7=1 


j= 

j^m 


and 


M 




{Zm, Zj)CjUe{Zj)€ F {Zm, Zj)K^{Zj)CjU^{Zj)Vol{^j). 


j=i 

j^m 


j^m 


Let US estimate the quantities A and B. 

2.5.1 Estimate of A 

By the decomposition of fl, 11 := \ we have 


(2.37) 


(2.38) 


r ’ r 

^‘^{zm,y)Ke{y)C^{y)U^{y)dy='^ {z^,y)K^{y)C,^{y)U^{y)dy. (2.39) 

J Vt j 1 J Q; 


Hence, A := f {zm,y)K^(y)Ceiy)Ueiy)dy 


E 

j=i 

j^m 


{Zm, Zj)Keizj)CjUeizj)Vol{nj) - / {z^, y)Keiy)Ceiy)Ueiy)dy 


(2.40) 


For 1 ^ TO, we have 



^‘^{zm, y)K^{y)C^{y)U^{y)dy 


r“(z™, zi)K,{zi)CiU,{zi)Vomi) 
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= K,{zi)Ci[ [T'^{zm,y)U,{y)-r‘-{z^,zi)U,(zi)]dy. (2.41) 

JQ, 


We set fizm,y) = r‘^{zm,y)U^{y) then every component f^{zm,y) of f{zm,y) satisfies 

fi{Zm,y) - fi{Zm,Zl) = {y - Zl)Rl{Zm,y) 

where 

R\{zm,y) = f Vyf,{zm,y - liiy - zi))d/3 
Jo 

pi 3 


do 

»1 3 

/ U,^j{y - P{y - zi))dp 

do 

»1 3 

+ / ~ (diy - zi))[yyU,j{y - P{y - zi))] dp. 

Jo 


From (2.6l and from Section 2.2 we derive for / m 

c 


|r“j(2:m,y-/3(y-2:i))l < 


i 1 
e3 


and \VyT‘p^j{zm,y - Piy - zi))\ < 


Att 


47rn2 ( ^ 


where c depends only on w and some universal constants. Then 


(2.42) 


\Ri{z^,y)\ < -^(JT['\U,{y-P{y-Zi))\dp+r\WyU,{y-P{y-Zi))\dp). ( 2 . 43 ) 

ZTT ne 3 \ne^ Jo Jo J 


Then, for I ^ m, (2.41) and (2.43) and observing that Ci is a constant matrix in fl;, imply the estimate 




T‘^{zjn,y)K^{y)C^{y)U^{y)dy - r“(zm,z/) K^{zi)CiUe{zi)Vol{Q.i) 


< /n, [ fo - Idiy - zi))\ dp] \y-zi\dy 

TT n^e3 * L 


< [K,(zi)]Ci 4 < 

„ 2 j n 2 


iTT ne3 Jo, Jo 


|y- z/My 

(2.44) 


for a suitable constant ci. 

Regarding the integral r‘^{zm,y)C^{y)Ue{y)dy we do the following estimates: 


r‘^(2m, y)Ke{y)C^{y)Ue{y)dy 


2.33l ^2Re(Zm)C^m [T (^mj J/) Ioo^^?/ 

'^2Ke{Zm)Cm(^fB(zm,r) |z„-y| ^ 


< e 

l 2 . 6 fc 47r 


i\B{zm,r) \Zm 2/1 




(here, -r G L^{B{z^,r)),r < 

^ \zm - y\ 2 / 
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< 


= 27rr 


{zm)Cm(o-(S^ /o^ ^ ds+^Vol(nni\B{zjn,r))'j 


e - 7rr 

3 


£zC2CmKe (^m) 


— :lin{r,e) 

< i^C2K^{zm)Cmlm(r^,e), 


— ( 3 ^) ^ e3 . 


here is the value of r where lm{r, e) attains maximum. 

1 

drlmij, e) = 0 ^ 47 rr- - — -irr = 0 => rj, = 

lm{r,,e) = 2 ^ (| 7 r) * el + (;^) ^ el - (| 7 r) * ei 

(2.45) 


From (2.401, we can have 


|^|oo<| [ T‘^{zm,y)K^{y)Ce{y)Ueiy)dy\c 


■E 

i=i 

j^m 


\T‘^{z^, Zj)K^{zj)CjUe{zj)Vol{nj) - / r‘^{z,n,y)K^{y)C^{y)Ue{y)dy\c 


which we can estimate by 
[E-i] 

|>i|oo < ^ 2[(2n+l)3-(2n-l)3] 


\T'^{z^,z,)K,(zj)C,U,{Zj)Vol{nj) 


and then 
Finally 

2.5.2 Estimate of B 

M 


r“( 2 m, y)K^{y)C^{y)Ue{y)dy\c 

l^loo ^ ^3^d^max[^^ 4“ ^^]- 
1^1 ^ C4CKjjifix^^ • 


y)K^{y)C^{y)Ue{y)dy\^. 


^ Zj)CjC/e(Zj)e ^ r {Zm-, Z'j )d^e {z,j)Cj Ue {Zj)Vol{Vtj ) 


i=i 

j^m 


t=l 

[Kd^^)] [e-^][K,{z,)] 


j^m 


= X! Zl)ClUe{zi)e + ^ ^ T'^{Zm,Zl)ClUe{zi)e- ^ {Zm, Zj)K^{Zj)CjUe{Zj)Vol{Qj) 


1^1 

l^m 

ZlGQrn 


i=i ^=1 

j^m zi£Qj 


i=i 

j^m 


IK,M] h-b 

= Cme ^ T'^{z^,Zi)Ue{zi) -\- ^ C'je[( ^ {Zm, Zi)Ue{zi)) -T'^{Zm,Zj)[K^{Zj)]Ue{Zj)], 


/=1 

Iz/zm 

Zl^^m 


i=i /=i 

j^m ziGQj 
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since Vol{VLj) = and Ci = Cj, for I = 1,..., K^{zj). We write, 


and 


E{ := ^ T^{zm,zi)U,[zi) 


(2.46) 


1=1 

l^m 

ZiGClm 


[K’e {Zj )] 


^2 := [( E 


Z=1 
zi Gflj 


[KAzj]] 

E (r“(z„,Zi)t/e(z/) - T‘^{Zjn,Zj)UeiZj)). 


(2.47) 


/=! 

2 ^ 


We need to estimate Cm^El and CjcE^. 

j^m 


Let us write, f'{zm,y) ■= (.Zm,y)Ue{y)- For zi G ftj, j ^ to, using Taylor series, we can write 

f'{Zm,Zj) - f'{Zra,Zl) = [Zj - Zi)R! [Zm] Zj, Zi), 

with 


R'{Zm]Zj,Zi) = / Vyf'{Zm,Zj - I3{zj - Zl))dl3. 

JO 


(2.48) 


By doing the computations similar to the ones we have performed in ( 2.42|[2(43 ) and by using Lemma 2.5 
we obtain 


^ ^ ^ C 4 C 


(2.49) 


t=i 

j^m 


One can easily see that. 


471 d= ■ 


(2.50) 


2.5.3 End of the approximation by the algebraic system 


Substitutio n of (|2.39 l in (2.38) and using the estimates (2.44) and (2.45) associated to A and the estimates 
(2.49) and (2.50) associated to B gives us 


Ue(z^) + E r^(z^,Zy)CyU,(z,)e = UHzm,0) + O (ciK^axe^J + O (2.51) 

j=i V ^ / 

j^m 
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We rewrite the algebraic system (1.8) as 


M 


= U\z^) 


i=i 

j^m 


where we set ■= —CraQm, recalling that Cm = Cm £• 


Comparing this system with (1.8) and by using Lemma 2.3 we obtain the estimate 


(2.52) 


Taking the difference between (2.51) and (2.52) produces the algebraic system 

M 

{u^^m ~ Ue{Zm)) + ^ i^rm Zj)Cj(UeJ — Ue{Zj))e = O {CKmax{^^ + ■ 

i=i 

j^m 


M 


Y,{U,,m-U,{Zm)) = O {CKmaxM{e^ + . (2.53) 

m—1 

For the special case d = e‘, M = 0{e~^) with t > 0, we have the following approximation of the far-field 


from the Foldy-Lax asymptotic expansion ( |l.6[ ) and from the definitions Uc^m '■= —Cm^Qm and Cm ■= Cm^, 
for m = 1,..., M: 


M 


inclU^{x,e)-x = 

i=i 


X e 


(2.54) 


-\-0 -j- ^2 —5i+3ia _|_ ^3—9t-\-6ta _|_ ^l — 2ta _|_ ^1 —3ia _|_ ^2 —5i+2ta^ 


Consider the far-field of type: 


/e "‘=p''"^K^{y)C^{y)Ue{y)dy--x(^x) [ e " K^{y)C^{y)Ue{y)dy. 
47rc‘ Jq AirCg Jq 


corresponding to the scattering problem (2.34) and set 


1 


and 


Uc„p{^^^) ■=-^:^ix(gix) e ' K^{y)C^{y)Ue{y)dy 


u^.A£,o) :=-^{i-x^x) [ e-'^s--yK,{y)cMU3y)dy 


47rc? 


(2.55) 

(2.56) 


Taking the difference between (2.54) and (2.55) we have: 
47rc2([/~(5,0)-t/g),p(x,0))-x 

. M 

Jn 

-\-0 (c -|- ^3—9i+6iQ _|_ ^l — 2toi _j_ ^l—3ta _j_ ^2—5t+2ia\ 


= E 


[e-i] [K,iz,)] 


— l — x 

e ""p 


i=i 


K,{y)C^{y)U^{y) -xdy-^ ^ 


-'i - X ‘ Zl 7-7- /V 

e “p CiU^i ■ xa 


j=i 1=1 

ziGflj 
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-\-0 -(- g2—5i+3ia _j_ ^3—9t-\-6ta _|_ ^1 —2fa _|_ ^l—3ta _j_ ^2—dt-\-2ta'j 


= 

i=i 


' fij '- 


Ue{y) ■ X — e ‘=p '^Ue(zj)-x dy 


3 = 1 


[Ke {Zj )] 




^ -l — X-Zl tt 

X — e Ue 


1^1 
_ Zl ^ O j 


[K’g (zj )] 




Z =1 

2:z GHj 


_|_(^ _|_ ^2—5t+3ia _j_ ^3—9i+6ia _|_ ^1 —2ia _|_ ^l — 3ta. _|_ ^2—5£+2iQ;^ 


= E / 


j=i 


e “"p 


^U^{y)-x — e pp '^Ug{zj)-x 


dy 


[KA^j)] 


E E (e -i-e' 

Z^l 

•S'J ^ ^2 -i 


M 


'■"'C/,(zz) • i) + E e [C/,(z,) - t/,,,] • i 

i -1 


_|_(^ _|_ ^2—5t+3ia _j_ ^3—9t-\-6toi _|_ ^1 — 2^0; _|_ ^l — 3ta. _|_ ^2—5t+2iQ;^ 


||2.53J E ^<^(^3)^: 

j = l 


n, ‘- 


-i — x-y 

e ""p 


f7e(?/) • i - e ^P ' ’UeiZj) ■ X 


dy 


+ E^je E (® *-=p""'C/.(2/)-i)+0(c2i^^a.(e5+ei-‘)) 


i=i z=i 

zi^Qj 


(2.57) 


_|_(^ _|_ g2—5Z:+3ia _j_ ^3—9i+6iQ _|_ ^1 —2fa _j_ ^l—3ta _j_ ^2—5t+2ia^ 

Now, let us estimate the difference Kt{zj)Cj e~^'^^'^U^{y) — U^{zj) dy. Write, fi{y) = 

e *pp“ ^U^{y). Using Taylor series, we can write 

fi{y) - /i(u) = iy- Zj)Rj{y), 


with 


{Rj{y))k,i = f Vy{fi)k,i{y-P{y-Zj))dl3 

J9 


'*1 r 




- -I k 


dp 


"1 r 




Ue{y - P{y - zj))dp 




[V,t7,(y-/3(y-z,))], dp. 


(2.58) 


We have V„e ^ = —i—xe '“"p^ ^ then 


|7?j(y)|oo < (— [ \Ue{y - p{y - Zj))\oodP + [ \WyUPy - P{y - Zj))\^dp'\ . (2.59) 

\Cp Jo Jo / 


Using (2.59) we get the estimate 


i=i 


/n,- '- 


e "p '%y)Ue{y) - e "p ' ^U^{zj) dy 


< 


(2.60) 
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'^K,{zj)Cji^ \y- 


j=i 


Cp Ja. 


\Ueiy - Piy - Zj))\dPdy 


+ J2K,{zj)Cj \y-Zj\\^vUe{y-ld(y-Zj))\^dl3dy 


^ ^ ^ ^e{Zj')CjCi 66 ^ “t“C 5 j^ Kfnax^^l ^ ^ “ 1 “ C 5 j . 

j=l 


In the similar way, using (2.53), we have 


[e-i] [K,iz,)] 

j=i 1=1 

^ J 


< 


O {^KmaxCiei + e^-^)). (2.61) 


Using the estimates (2.60) and (2.61) in (2.57), we obtain 


A7rcl{U^{x,9)-U^^^p{x,e))-x = 


Since Vol{V,) is of order e ^diam(i?m) + |j , and d is of the order e*, we should have t > ^. 
Hence 

t>-; 1 — t>0; 3 — 9t + 6 ta >0; 1 — 3ta >0; 2 — 5t + 2ta > 0; 

O 

Hence for | < t < 1, we have 

ta < -] 2 — 5t + 2ta >0; 3 — 9t + 6 ta > 0. 

O 

Equating 2 — 5t + 2ta = 3 — 9t + 6 ta, we find that at = t — ^ and then 2 — 5t + 2ta = 3 — 9t + 6 to = | — 3t 
and 1 — 3af = j — 3t. In addition, since at < |, then t < To make sure that | — 3t is positive we should 
have t < 5 . Hence, the error is 


O ( Kjnax^^Cci (- 1 - C 5 ) ) + 0{C{C + l)CKmax{^^ + 6 ^ *)) 


-\-0 ^2 —5i+3ta _|_ ^3—9t-\-6ta _|_ ^l — 2ta _|_ ^l — 3ta. _|_ ^2 — 5t-\-2ta\ 


o 


_|_ — i _j_ ^3—9t-\-6a.t _|_ ^l — 3<y.t _|_ ^2—5i+2ia\ 


(2.62) 


47rcl[U^{x,9)-U^^Jx,e)]-x = O ^ ^ 


(2.63) 


Doing similar computations for the s-parts of the far fields as it was done for the p-parts, between (2.54- 


2.62), we obtain the following estimates 

47rcl[Urix,0)-U^^^M0)]-x^ = o(el+ei-3y, l<t 


1 

^ 2' 


(2.64) 
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2.6 End of the proof of Theorem |1.2 

Combining the estimates (|2.63[) and (2.37), we deduce that 


iiTcl [c/“(i,d) - U^^ix,B)] ■ X = 0(e-M7. 1-3*}), e « 1, 


uniformly in terms of x,0 € S^. 


Similarly, by combining the estimates (2.64) and (2.37), we deduce that 


Ancl [t/r (i, 0) - U^^x, d)] • = 0(. 


-L _ i. f-3*} 


uniformly in terms of x, 9 G S^. 


1 1 

7, <t< 

3 “ 2 


), e « 1 , o<t<7; 


(2.65) 


( 2 . 66 ) 


3 Justification of Corollary 1.3 


For an obstacle of radius e, 5((())(s) := and 'D{(j)){s) := Similarly, 

we set 5 g(()i)(s) := jQjj^Gp{s,t)4>(t)dt and 'Dg{4'){s) := 4>(t)dt. We see that VF«,(a;, z) := 

Gp{x, z) — r“(a:, z) satisfies 

{^ + uj^)W^=u:^{l-p)T^, inM^ (3.I) 

with the Kupradze radiation conditions. Since r‘^(-,z), z € is bounded in LP{fl), for p < 3, by interior 
estimates, we deduce that bF(-,z), z G is bounded in W^’P(n), for p < 3, and hence, in particular, the 
normal traces are bounded in L^{dD^). Then we can show that the norms of the operators 


and 


are of the order 0(e) at least. 


5g - 5 : {L^{dD,)f ^ {H\dD,)f 


Vg-V-. {L^{dD,)f ^ {L^{dD,)f 


(3.2) 


(3.3) 


1. Using these properties and arguing as in [Challa fc Sini, 2015 , we derive the asymptotic expansions 
(1.25)-(1.26). Indeed, apart from the computations done in [Challa &: Sini, 2015| , the main arguments 
needed to extend those results to the case of variable density is the Fredholm alternative for the 
corresponding integral operators and the application of the Neumann series expansions. After splitting 


Gp as Op = r“ + (Op — F"^), these two arguments are applicable as soon as we have (3.2)-(3.3). 


2. The justification of the invertibility of the algebraic system (1.8) depends only on (1) the distribution 


of small bodies and (2) the background medium through the singularities of the fundamental solution 
(of the form jr“(s,t)| < c|s — t|“^). However, this type of singularity is true for general background 
elastic media|j Then the same arguments can be used to justify the invertibiliy of the algebraic system 
(1.27). Using the above mentioned decomposition of the Green’s function Gp, the properties of the 


Lippmann Schwinger integral equation are also valid replacing r“ by Gp, and hence the results in 
section 2.4 are valid. Finally, and again using the decomposition of Gp, the computations in section 


2.5 can be carried out using Gp. 


4 The elastic capacitance 

We start with the following lemma on the symmetry structure of the elastic capacitance. 


®Of course, it can be justify using the decomposition Gp = r“ + (Gp — r“) with the singularity of r“ and the smoothness 
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Lemma 4.1. Let C := ijgjy be the elastic capacitance of a bounded and Lischitz regular set 

D and C* be its adjoint. Then 

C = C*. (4.1) 


Proof of Lemma j.l We know that the matrix a := solves the invertible integral equation 


IdD^o{s,t)a{t)dt = I 3 , or precisely JgjjTQ{s,t)ai{t)dt = e^, where at := and Ci is the column 

of I 3 . Let a be any constant vector in then the vector a a satisfies 



ro(s,t) {cr{t) a) dt = a. 


We set := JgjjTois,t) {cr(t) a) dt. Then ip°' satisfies the problem = 0, in and = a on dD. In 

addition, we have the jump relation = a{t) a on dD where d^u := A(V • M);/ + /i(Vu + 

is the elastic conormal derivative. Hence 


dy+‘p°' — d^-ip^ dt = 


aft) adt = C a. 


JdD JdD 

Now, let a and b be arbitrary constant vectors in To both a and b, we correspond and as above. 
Using the Green formulas inside and outside of U, we deduce that 


(C a,b) = f (i9i/+iy5“ - d^-ip°-(t)) ■ ip^{t) dt = f {d„+(p^ - dy-ip^) ■ ip°‘dt = {C b, a) = (a, C b) 
JdD JdD 

recalling that every quantity here is real valued. 


The next lemma describes the elastic capacitance of a given bounded and Lipschitz regular domain with 
the one of its image by a unitary transform. 

Lemma 4.2. Let TZ = (rim) be a unitary transform in D be bounded Lipschitz domain in d = 2, 3 
and D = TZ{D). Let C and C be the corresponding elastic capacitance matrices due to the density matrices 


a and a, as defined in (1.9), respectively. Then we have C = TZ C TZ 


-1 


Proof of Lemma j.2 First recall the relation T® o 7?.(^, rf) = 7^r°(^, rfjTZ see [Ammari &: Kang, 2004[ 
Lemma 6.11]. From (1.9), we have that 



Sdb r°(l,d)CT(i7)dd 

= I, i) G au 


IdD (r° 0 (C, d) (a- 0 TZ) fq)d'q 

= I, 7? G au 


fdD ^r°(C, ri)TZ-^d 0 TZ{r])dr] 

= I, 7 ? G au 


IdD (cr 0 7^) (?7)d?7 

= 77 G an 


Jg^r°{f,rj)TZ-^{doTZ) {r^)TZdig 

= I, 7 ? G au. 

Now from the uniqueness of solutions of (|1.9|), we deduce that 

7^-l (a•o7^) f)TZ = 


(CTo7^) (•) = 7^cr(•)7^■^ 

From the definition of the capacitance, see (11.91), and (14.3|), we have 


C= a(fj)dfj = / {d o TZ) fq)d'q = / TZafq)TZ ^dri = TZCTZ . 
JdD JdD JdD 


(4.2) 

(4.3) 

(4.4) 
□ 


Proposition 4.3. Let TZ = (rim) be a unitary transform in d = 2,3, and let D be a bounded Lipschitz 
domain in d = 2,3, and D = TZ{D). Let C and C he the corresponding elastic capacitance matrices due 
to the density matrices a and a, as defined in (1.9), respectively. 
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1. 2D-case. If the shape of D is rotationally invariant for any rotation hy one angle 0 0 ,tt, then C is 

a scalar multiplied by the identity matrix. 

2. ZD-case. If the shape of D is rotationally invariant for any two of the rotations around the x,y or 
z axis by one angle 0 7 ^ 0, tt and a 7 ^ 0, tt respectively, then C is a scalar multiplied by the identity 
matrix. 

Proof of Proposition |//.3| In 2D case, the rotation matrix by an angle 9 is given by 

= 


cos 9 — sin 9 
sin 9 cos 9 


(4.5) 


As the shape is invariant by this rotation then C = C. Since IZ is unitary then TZ ^ = IZJ and then (4.4) 
implies 


C 


= C = TZCTZ^ 
cos 9 — sin 9 
sin 9 cos 9 


Cii <^12 \ ( cos 9 sin0 
C 21 C 22 ) \ — sin0 COS0 

COS0 — sin0 \ / Cii COS0 — C'i 2 sind Cn sin 0 + (712 cos 0 
sin 9 cos 9 J y C 21 cos 9 — C 22 sin 9 C 21 sin 9 + C 22 cos 9 

Cii cos^ 9 + C22 sin^ 9 — {C12 + C21) sin 6 * cos 6 * C12 cos^ 9 — C21 sin^ 9 + (Cn — C22) sin 0 cos 0 

C21 cos^ 9 — C12 sin^ 9 + {Cn — C22) sin 0 cos 0 Cn sin^ 9 + C22 cos^ 9 + {C12 + C21) sin 0 cos 0 


(4.6) 


We deduce from (4.6) and the symmetry of matrix C the following relations: 

Cii = Cii cos^ 9 + C22 sin^ 9 — 2 Ci 2 sin 9 cos 9 

C22 = Cii sin^ 9 + C22 cos^ 9 + 2Ci2 sin 9 cos 9 
C12 = C\2 cos^ 9 — C12 sin^ 9 + (Cn — C22) sin 9 cos 9 . 


We rewrite (4.7) and (4.9) respectively as 

{Cii — C 22 ) sin^ 9 + 2 Ci 2 sin cos 9 — Q 

{Cii — C 22 ) cos 0 sin 0 — 2 Ci 2 sin^ 0 = 0 
Taking sind as a multiplicative factor in ( 4.10| ) and (4.11) we see that if 0 7 ^ 0, tt, i.e. sin0 7 ^ 0, then we have 

Cii = C 22 and C 12 = C 21 = 0. (4-12) 


(4.7) 

(4.8) 

(4.9) 

(4.10) 

(4.11) 


Let us now consider the 3D case. First, let us assume that the shape is invariant under the rotation 
about the x — axis and an angle 0 7 ^ 0, tt. This rotation matrix is given by 


1 0 0 

7^ = I 0 cos 0 — sin 0 

0 sin 0 cos 0 


Since TZ is unitary, TZ ^ = TZ^ , then (4.4) gives us; 
C = C = TZCTZ^ 


1 

0 


/ Cn 

C12 

Ci3 

0 

COS0 

— sin 0 

C21 

C22 

C23 

0 

sin0 

COS0 / 

\ C31 

C32 

C33 



(4.13) 


0 0 
cos 0 sin 0 

— sin 0 cos 0 
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(4.14) 


where 


Til 

Ti2 

Ti3 

T21 

T 22 

T 23 

T31 

T32 

T33 


1 

0 

0 \ f Cn 

(7i 2 cos 9 — 

(7i 3 sin 6* 

(7i 2 sin 9 + (7i3 cos 9 

0 

cos 9 

-sin 9 C21 

C22 cos 9 - 

C23 sin 9 

C22 sin 9 + C23 cos 9 

0 

sin 6* 

cos 9 J \ C31 

C32 cos 9 - 

C33 sin 9 

C32 sin 9 + C33 cos 9 

Til 

Ti2 

Ti3 \ 




T21 

T22 

T23 , 




T31 

T32 

T33 / 





Cii; 

(712 cos 0 — (7 i 3 sin 0; 

(7i 2 sin 9 + (7i3 cos 0; 

(721 cos 0 — (731 sin 0; 

(722 cos^ 9 + C33 sin^ 9 — {C22 + C33) sin 9 cos 9 ; 

C23 cos^ 9 - C32 sin^ 9 + {C22 - C33) sin 9 cos 9 ; 

C21 sin 9 + C31 cos 9 ] 

C32 cos^ 9 - C23 sin^ 9 + {C22 - C33) sin 9 cos 9 ; 

C22 sin^ 9 + C33 cos^ 9 + ((723 + C32) sin 9 cos 9 ; 


(4.15) 


We observe the equality of the 2x2 matrices: 


C22 

(723 

\ - 1 

^ T 22 

T 23 \ 

C32 

CO 

1 - 1 

y T 32 

T 33 J 


(4.16) 

C22 cos^ 9 + C33 sin^ 9 — {C22 + C33) sin 9 cos 9 C23 cos^ 9 — C32 sin^ 9 + {C22 — C33) sin 9 cos 9 

C32 cos^ 9 — C23 sin^ 9 + {C22 — C33) sin 9 cos 9 C22 sin^ 9 + C33 cos^ 9 + {C23 + C32) sin 9 cos 9 

These are similar to the matrices we obtained in the 2D case. Hence we deduce, as in the 2D case, that 

C22 — C33 and C23 = C32 = 0. (4.17) 

To show that (7 is scalar multiplied by the identity matrix we need to prove that Cn = C22, for instance, 
and Ci3 = C31 = 0 . For this purpose, we use another rotation. Taking the rotation around the z-axis Qby 
one angle a 7^ 0, tt and proceeding as we did for the rotation about the x-axis, we show that 


(7ii — C33 and (7i3 — (73i — 0. 


13 — 0^31 


(4.18) 

□ 


From the above analysis, we have the following remark: 

Remark 4.4. 1. For the spherical shapes, in particular, the capacitance is a scalar multiplied by the 

identity matrix. 

2. Ellipsoidal shapes are invariant only under rotations with angle -k (or trivially Q). For these shapes, 
the capacitance might not he a scalar multiplied by the identity matrix but a diagonal matrix instead. 
To justify this property, the arguments in \Ammari & Kang, 2004^ can be useful. 

numbib] imamat 

^ We can use also the rotation around the y axis. 
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